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Abstract — We consider the problem of network coding across 
three unicast sessions over a directed acyclic graph, when each 
session has min-cut one. Previous work by Das et al. adapted 
a precoding-based interference alignment technique, originally 
developed for the wireless interference channel, specifically to 
this problem. We refer to this approach as precoding-based 
network alignment (PBNA). Similar to the wireless setting, PBNA 
asymptotically achieves half the minimum cut; different from the 
wireless setting, its feasibility depends on the graph structure. 
Das et al. provided a set of feasibility conditions for PBNA 
with respect to a particular precoding matrix. However, the set 
consisted of an infinite number of conditions, which is impossible 
to check in practice. Furthermore, the conditions were purely 
algebraic, without interpretation with regards to the graph 
structure. In this paper, we first prove that the set of conditions 
provided by Das. et al are also necessary for the feasibility of 
PBNA with respect to any precoding matrix. Then, using two 
graph-related properties and a degree-counting technique, we 
reduce the set to just four conditions. This reduction enables 
an efficient algorithm for checking the feasibility of PBNA on a 
given graph. 

I. Introduction 

Network coding was originally introduced to maximize the 
rate of a single multicast session over a network |TJ Q 
||3). However, network coding across different sessions, which 
includes multiple unicasts as a special case, is a well-known 
open problem. For example, finding linear network codes for 
multiple unicasts is NP-hard |4|. Thus, suboptimal, heuristic 
approaches, such as linear programming |5| and evolutionary 
approaches (6|, are typically used. Moreover, while it has been 
shown that scalar or vector linear network codes might be 
insufficient to achieve the optimal rate [7], only approximation 
methods [8] exist to characterize the rate region for this setting. 

In this paper, we consider the simplest inter-session linear 
network coding scenario: three unicast sessions over a directed 
acyclic graph, each session with minimum cut one. Das 
et al. |9j applied a precoding-based interference alignment 
technique, originally developed by Cadambe and Jafar fT0[ 
for wireless interference channel, to this problem; we refer to 
this technique as precoding-based network alignment (PBNA). 
In a nutshell, PBNA (i) simulates a wireless channel through 
random network coding p) in the middle of the network 
and (ii) applies interference alignment at the edge, i.e., via 
precoding at the sources and decoding at the receivers. This 
way, it greatly simplifies the network code design, while it 
guarantees that each unicast session asymptotically achieves a 
rate equal to half of its minimum cut (9). 

An important difference from the wireless interference chan- 
nel is that, in our problem, there may be dependencies between 



elements of the transfer matrix introduced by the graph struc- 
ture, which make PBNA infeasible in some networks (TTJ. As 
a first step, Das et al. [9 1 provided a set of feasibility conditions 
for PBNA, and proved they are sufficient for the feasibility 
of PBNA with respect to a particular precoding matrix. One 
important limitation is that the set consists of an infinite 
number of conditions, which makes it impossible to check 
in practice. Another limitation is the lack of consideration 
of graph structure, which turns out to be the reason for the 
significant redundancy in the set of conditions. Ramakrishnan 
et al. [11 1 conjectured that the infinite set of conditions can be 
reduced to just two conditions. Han et al. fl2) proved that the 
conjecture holds for three symbol extensions; however, this 
result cannot be generalized beyond three symbol extensions. 

In this paper, we make the following contributions. First, 
we prove that the set of conditions provided in [9| are also 
necessary for the feasibility of PBNA with respect to any valid 
precoding matrix. Then, by using a simple degree-counting 
technique and two graph-related properties, we greatly reduce 
the set to just four conditions; two of them turn out to have an 
intuitive interpretation in terms of graph structure. Finally, we 
present an efficient algorithm for checking the four conditions. 

The rest of this paper is organized as follows. In Section 
II, we present the problem formulation. In Section III, we 
summarize our main results. In Section IV, we discuss the 
graph-related properties that are key to the simplification of 
the conditions. In Section V, we prove and discuss our main 
results regarding the feasibility condition of PBNA. In Section 
VI, we present an algorithm for checking the condition. In 
Section VII, we conclude the paper. The Appendices provide 
details on the proofs that were outlined or omitted from the 
main part of the paper. 

II. Problem Formulation 

The network is a delay-free directed acyclic graph, denoted 
by G — (V, E), where V is the set of nodes and E the set of 
edges. Without loss of generality, each edge has capacity one, 
i.e., can transmit one symbol of finite field F 2 ™ in a unit time. 
For the zth unicast session (i € {1, 2, 3}), let Si and di be its 
sender and receiver respectively, and Ri its transmission rate. 
Every edge e€£ represents an error free channel. We assume 
that the minimum cut between s,; and di is one. Let Xi be the 
source symbol transmitted at Sj and Zi be the symbol received 
at di. We further extend G as follows: For the ith unicast 
session (i £ {1, 2, 3}), we add a virtual sender s[ and a virtual 
receiver dl i and two edges Oi — (s^Si) and Tj = (d,,^). The 



extended graph is denoted by G — (V',E r ). For e G E', let 
head(e) and tail(e) denote its head and tail respectively. 

In the middle of the network, we employ random network 
coding [3|. The symbol transmitted along e e E', denoted by 
Y e , is a linear combination of incoming symbols at tail(e). 



Y, = 



Xi Ife = (Tj; 

K E h ead(e>)=taU(e) x e> e Y e > Otherwise, 
where variable, which takes values from F 2 >> 



and 



represents the coding coefficient used to combine the incoming 
symbol along e' into the symbol along e. We group all coding 
coefficients x e / e 's into a vector x, called the coding vector of 
G' . The network acts as a linear system: the output at d\ is 
a mixture of source symbols, 2, = Y^j=i m ij{ K )-^j' where 
rriij (x) £ F 2 ™ [x] is the transfer function from s'j to and 
can be written as follows (2): 



my(x) 



where 'Py is the set of paths from s'j to c^, and t(P) is the 
product of coding coefficients along path P. We assume that 
all my(x)'s are non-zeros, which is the most challenging case. 
Indeed, as shown in Section [V] when some my (x) (i 7^ j) 
is zero, the feasibility condition of PBNA is significantly 
simplified due to reduced number of interferences. 

At the edge of the network, we apply interference alignment 
1(9) (lOj via precoding at senders and decoding at receivers. 
Let Xi = (Xl, • • • , denote the input vector at sender 

s'i, where fej = n + 1 for i = 1 and fcj = n otherwise. We 
use precoding matrix Vj to encode Xj into 2n + 1 symbols, 
which are then transmitted via 2n + 1 uses of the network 
(time slots). The output vector at d! i is 

v^ 3 



M,,V,X, 



where My- is a (2n + 1) x (2n + 1) diagonal matrix with the 
(k, k) element being my(x fe ), where x fc represents the coding 
vector for the fcth use of the network. Vi is a (2n+l) x (n+l) 
matrix, and V2, V3 are both (2n+l) x n matrices. Vj can still 
contain indeterminate variables. Let £ denote the vector of all 
variables in x 1 , • • • ,x 2n+1 and V 1; V 2 , V 3 . PBNA requires 
that the following conditions are satisfied for some values of 

span(Mi 2 V 2 ) = span(Mi 3 V 3 ) 
span(M 23 V 3 ) C span(M 2 iVi) 
span(M 32 V 2 ) C span(M 3 iVi) 
rank(M 11 V 1 M 12 V 2 ) = 2n + 1 
rank(M 2 iVi M 22 V 2 ) = 2n + 1 
rank(M 3 iVi M33V3) = 2n + 1 

Condition s^i guarantees that all the interferences at d! i 
are aligned, i.e., mapped into the same linear space; while 
condition SSi ensures that all source symbols for the ith 
unicast session can be decoded. These conditions ensure 
that we can achieve a rate tuple (Ri, R2, R3) = R* t = 



In this case, we say that R* is AfA-/eas/bZe[|] 

Previous work |9| |11| [12] only considered the feasibility 
of PBNA under a particular precoding matrix, i.e., in 
Eq. (|3j, which was first introduced in | [T0| . To address this 
limitation and characterize the feasibility of PBNA for any 
precoding matrix, we reformulate srf\, 3/2 , ^3 and SS\ , S82 , ^3 
without any assumption about the structure of precoding 
matrix. First, we reformulate as: 

si[ : M 12 V 2 = M 13 V 3 A 
s% : M 23 V 3 = M 21 V!B 
44 : M 32 V 2 = M 31 V!C 

where A is an n x n invertible matrix, and B and C are both 
(n + l) x n matrices with rank n. ^/ 2 , s^l can be further 
condensed into a single condition: 



TViC = ViBA 



(1) 



where T = Mi 2 M 2 " 1 1 M 23 M 32 1 M 3 iM ] ; 3 1 . Finally, conditions 



l 32 

are reformulated as: 

lMf) = det(Vi 
MO=det(V 1 
^3(0 = det(Vi 



P1V1C) ^0 
P 2 V!C) ^0 
P^xCA- 1 ) ^0 



M 31 M 2 ' 1 1 M 22 M 



32 



where Pi 

and P 3 = M 12 M 3 " 2 1 M 33 M 13 1 , and MO, MO, MO 
are rational functions in the field F 2 m(£). Define ip(0 = 
lli=iV'i(C)- We assume that F 2 ™ is sufficiently large such 
that if ip(0 is a non-zero rational function, there are values 
to £, denoted by £0, such that ip(£,o) 7^ 0. 

We also define the following rational functions: 



Pi(x) 



m 3 i(x)mi 2 (x) 



P8(x) = 



m 3 i(x)m 2 2(x) 



(2) 



m 1 i(x)m 82 (x) m2i(x)m 32 (x) 
mi 2 (x)m 3 3(x) , . _ m 3 i (x)m 12 (x)m 23 (x) 
m 32 (x)mi 3 (x) m 2 i(x)m 32 (x)mi 3 (x) 

Clearly, pi(x) and r)(x) form the elements along the diag- 
onals of Pi and T respectively. Hence, the following lemma 
holds: 

Lemma 1: R* is NA-feasible if and only if 1) Eq. ([T]) is 
satisfied, and 2) 3B' X , are satisfied. 

The fundamental design problem in PBNA is to find Vi 
such that all the conditions in Lemma [TJ are satisfied. Indeed, 
the major restriction on Vi comes from Eq. ([TJ. As shown 
in |9j, the construction of V\ depends on whether 77 (x) is 
constant. When r\ (x) is constant, and thus T = I 2 „+i, we set 
C = BA. Therefore, any arbitrary Vi can satisfy Eq. ([JJi. In 
contrast, when ij (x) is not constant, we can no longer choose 
Vi freely. flO) proposed the following precoding matrix, 
which has also been used by most of recent work (5J 1 1 1 1 

& 



V* = (w Tw ■ • • T"w) 



(3) 



I n+l 
V 2n+l 



+ 1 ' 2ra+l ' 2n+l ■ 



which approaches (|, 4, i) as n 



where w is a column vector of 2n + 1 ones. Meanwhile, we 
set A = I„, C consists of the left n columns of I n +i, and B 

'In this paper, we first consider the feasibility conditions of PBNA for a 
fixed value of n. Then, in the Main Theorem, we prove that the feasibility 
conditions of PBNA are actually irrelevant to n for n > 1 . 




(a) pi(x) = p 2 (x) = p 3 (x) 

7J(X) = 1 



Fig. 1: Examples of graphs where PBNA is infeasible 

the right n columns of I„+i; this construction satisfies Eq. ([lj. 
Note that the form of Vi is determined by A, B and C. With 
different A.B and C, we can derive different Vi; therefore 
the choice of Vi is not limited to just V^. 

As observed in 1 1 1 1, graphs can introduce dependence 
between transfer functions^ so that PBNA may be infeasi- 
ble. This is a fundamental difference compared to wireless 
interference channel, where channel gains can change inde- 
pendently and interference alignment is always feasible. Fig. 
[T] depicts some examples of graphs where PBNA is infeasible. 
In Fig. [jja), p;(x) = ?/(x) = 1 for i G {1,2,3}, thus 
P». = l2n+l> implying 33' 2 , & 3 are all violated. In Fig. 



which also violates This example 



snows that the conjecture proposed by Ramakrishnan et al. 
[jTTJ doesn't hold beyond three symbol extensions. Moreover, 
by exchanging s\ O s 2 and di ■<-» d 2 , we obtain another 
counter example, where P2(x) = 1 + »j(x), violating 3$ 2 , 

As a first step, |9| proposed the following set of conditions 
for PBNA^ For i e {1,2, 3}, 



■" (x,s *{iwi :/w ' 9( 



=)eF 2 .,W, 9 ( s )^0 (4) 



In (9), it was proved that if Eq. Q is satisfied, we can use 
to asymptotically achieve half rate through infinite number 
of uses of the network. Unfortunately, since Eq. (|4| contains 
an infinite number of conditions, it is impractical to verify. 
Moreover, since only one particular matrix was considered in 
j9), Eq. Q was only shown to be sufficient for PBNA. 

III. Overview of Main Results 

We now state our main results; proofs are deferred to 
Section V and to the appendices. Since the construction of 
Vi depends on whether r\ (x) is constant, we distinguish two 
cases. 

A. 7y(x) Is Constant 

In this case, we can choose Vi freely, and thus the feasi- 
bility condition of PBNA can be significantly simplified, as 
stated in the following theorem: 

Theorem 1: Assume r/(x) is constant. R* is NA-feasible if 
and only if p»(x) is not constant for each i G {1,2,3}. 

2 Dependence here means that one transfer function (namely Wii(x), 
corresponding to signal for the ith unicast flow) can be written as a rational 
function of other transfer (interference) functions. The exact functional form 
is dictated by Eq. |5l or Eq. 161. 

3 There is actually a smalldifference between Eq. |4j and the original 
formulation in (91, in which pi(x) is replaced by l/pi(x). It is easy to 
see that the two are equivalent. 



B. rj{~x) Is Not Constant 

In this case, we cannot choose Vi freely. Using similar 
technique as in [9|, we can rewrite Eq. Q as follows: Q 

ftW iS n = {^jy : /(*),$(*) G F 2m [z], 

f(z)g(z)^0,gcd(f(z),g(z)) = l, (5) 
d f <n,d g <n-l\ Vie {1,2,3} 



Note that, in contrast to Eq. (|4]), the above set of conditions 
guarantee that R* is NA-feasible for a fixed value of n. 

Next, we show that Eq. <|3j is also necessary for the 
feasibility of PBNA with respect to any Vi satisfying the 
conditions of Lemma [T] 

Theorem 2: Assume 77 (x) is not constant. R* is NA- 
feasible if and only if for each i G {1,2,3}, Pi(x) ^ S n . 

Finally, we greatly reduce S n to just four rational functions: 

Theorem 3 (The Main Theorem): Assume rj(x.) is not con- 
stant. For n > 1, R* is NA-feasible if and only if for each 
i G {1,2,3}, the following condition is satisfied: 



Pi(x)gS'= l,77(x),r,(x) + l, 



??(x) + 1 



(6) 



As shown in the Main Theorem, the feasibility conditions 
for PBNA are irrelevant to n for n > 1. This indicates that 
if PBNA is feasible for n = 2, then it is feasible for any 
arbitrary n > 1, and thus we can use PBNA to achieve half 
rate asymptotically. Otherwise, if PBNA is not feasible for 
some n > 1, PBNA doesn't even allow us to achieve any rate 
greater than ( | , 1 , | ) . 

The basic idea behind the Main Theorem is that we can 
compare the degree of a variable in Pi(x) with that of a 
rational function in S n . This technique enables us to reduce 
S n to the form {fj^pf^^y}- Thus, we only need to consider a 
finite number of rational functions, namely S' in Eq. (|6|. This 
enables an efficient algorithm for checking the feasibility of 
PBNA. The key for enabling this reduction lies in two graph- 
related properties, which we refer to as Linearization Property 
and Square-Term Property, as described in the next section. 

IV. Graph-Related Properties 

Our key intuition is that Pi(x) is not an arbitrary function 
but depends on transfer functions, as specified in Eq. pi. 
Therefore, Pi(x) has special algebraic properties, which can 
be exploited to simplify Eq. ^3). 

First note that all pj(x)'s are of the following general form: 

m Q f,(x)m P9 (x) 

^( X ) = H f V a,b,p,qe {l,2,3},a^p,b^q 

Furthermore, each path pair in Tab x V vq contributes a term 
in m„j(x)m M (x), and each path pair in V a q x V v b contributes 

4 Notation: For two polynomials f(x) and g(x), let gcd(/(a;), g(x)) denote 
their greatest common divisor, and df the degree of f(x). 




(a) o(e 2 ) > o(e 3 ) 
and o(ei) < o(e^) 



(b) o{e%) > o(ez) 
and o(ei) > oiea,) 



(c) o(e2) < o(e3) 



Fig. 2: The construction of H (in the proof of the Linearization 
Property) enabled by Lemma [2] (Pi is disjoint with P2) 



a term in m aq (x)m p b(x): 
m ab (x)m pq (x) = 



t(Pi)t(P 2 



m aq {x)m pb (x) = V t{P 3 )t(P 4 ) 

A. Linearization Property 

First, consider the following lemma, which provides an easy 
way to check whether pi(x) ^ {1,t?(x)} (as in Section VI). 
The intuition is that we can multicast two symbols from s' b , s' q 
to d' a , d' p by network coding if and only if the minimum cut 
separating s' b , s' q from d' a , d' p is greater than one J2|. 

Lemma 2: m ab (x)m pq (x) ^ m aq (x)m pb (x) if and only if 
there is disjoint path pair (Pi,P 2 ) G V ab X P p(? or (P 3 , P 4 ) G 

Pag ^ Ppb- 

Proof: See Appendix |A| ■ 
The first graph-related property states that Pi(x) can be 
transformed into its simplest non-trivial form (i.e., a linear 
function or the inverse of a linear function). The key to Lemma 
[3] is to find a subgraph H and consider h(x) restricted to 

h, i.e., h(* H ) = ::X H X;t H H y where x * re P resents the 

coding vector of H. Due to the graph structure induced by 
Lemma [2j we can always find H such that some variable 
x ee i appears exclusively in the numerator or the denominator 
of /i(x#). Thus, by assigning values to xh other than x ee i, 
we can transform /i(x#) into a linear function or the inverse 
of a linear function in terms of x ee i. Since h(xn) can be 
acquired through a partial assignment to x, this transformation 
also holds for the complete graph G. 

Lemma 3 (Linearization Property): Let h(x) = 
m „ t (x)m M (x) = 2jx) h h gcd ( u (x)Xx)) = 1. Assume 

m a ,(x)m pb (x) u(x) & V V n \ I) 

h(x) is not constant. Then, for sufficiently large to, we can 
assign values to x other than a variable x ee > such that u(x) 
and v(x) are transformed into either u(x ee /) = cix ee ' + c , 
v(x ee >) = c 2 or u(x ee ') = c 2 ,v(x ee >) = cix ee > + c , where 
Co, ci, C2 are constants in F 2 ™, and c\c 2 7^ 0. 

Proof: In this proof, given a path P and e, e' G P, let 
P[e : e'] denote the path segment along P between e and 
e', including e, e'. We arrange the edges of G' in topological 
order, and for e G E' , let o(e) denote e's position in this 
ordering. Moreover, denote /ii(x) ~ m ab (x)m pq (x), /i2(x) = 
m aq (x)m pb (x) and d(x) = gcd(/ii(x), /i 2 (x)). Let si(x) = 
^ and s 2 (x) = ^g. Hence gcd( Sl (x), s 2 (x)) = 1. It 




Fig. 3: Illustration of Square-Term Property. A term with x 2 ee , 
introduced by (Pi, P2) in the numerator of h(x) equals another term 
introduced by (P3, P4) in the denominator of h(x). 

follows u(x) = csi(x) v(x) = cs 2 (x), where c is a non- 
zero constant in F 2 ™ . By Lemma [2] there exists disjoint path 
pair (Pi, P a ) G P a6 X P p? or (P 3 , P 4 ) G P QO x P pfc . Now we 
consider the first case. 

Let (P 3 , P4) G V aq x P pb . Since Pi, P^ both originate at a b , 
and P 2 , P4 both terminate at t p , there exist ei G Pi n P4 and 
e 2 G P 2 n P4 such that the path segment along P4 between e\ 
and e 2 is disjoint with P1UP2. Similarly, there exist e 3 G P 2 H 
P3 and e4 G PiHPg such that the path segment between e 3 and 
e4 along P3 is disjoint with P x U P 2 . Construct the following 
two paths: P'l = P^at, : e x ] U P^d : e 2 ] U P 2 [e 2 : r p ] and 
PJ' = P 2 [a q : e 3 ]UP^[e 3 : e 4 ]UPi[e 4 : r a ] (see Fig. [2}. Let H 
denote the subgraph of G" induced by Pi U P 2 U P 3 ' U P" , and 
x// the coding vector of H. We will prove that the theorem 
holds for H. Note that since /ii(x#) and h 2 (xu) are both 
non-zeros, d(xn) 7^ 0. 

If o(e 2 ) > o(e 3 ) (Fig. [2j»-(b)), the variables in t(P 2 [e 3 : 
e 2 }) are absent in h 2 (xn)- We then arbitrarily select a 
variable x ee i from i(P2[e 3 : 62]), and write h\{xn) as 
f(x' H )x ee > +g(x' H ), where x' H includes all the variables in xh 
other than x ee >, and f(x' H ), g(x' H ) G F 2 m[x' H ]. Meanwhile, 
h 2 {xu) can be written as h 2 (x' H ) G F 2 ™[x^]. Clearly, 
x ee ' will not show up in d(xn) and thus it can also be 
written as d(x' H ) G F2m[x^]. We then find values for x^, 
denoted by r, such that f(r)h 2 (r)d(r) 7^ 0. Finally, denote 
c = cy(r)d _1 (r), c\ = c/(r)d _1 (r) and c 2 = c/i2(r)d _1 (r) 
and the theorem holds. 

On the other hand, if o(e 2 ) < o(e 3 ) (see Fig. |2|c)), the 
variables in i(Pi[ei : e 4 ]) are absent in h 2 (xn)- We then 
select a variable x ee i from t(Pi[ei : 64]). Similar to above, 
it's easy to see that u(x) and u(x) can be transformed into 
c\x ee i + Co and C2 respectively. 

For the case where there exists disjoint path pair (P 3 , P4) G 
VagXTpt,, we can show that u(x) and i>(x) can be transformed 
into C2 and cix ee ' + Co respectively. ■ 

B. Square-Term Property 

The second graph-related property is stated in Lemma[4j the 
coefficient of x\ e t in the numerator of h(x) equals its counter- 
part in the denominator of h(x). Thus, if x 2 ee , appears in the 
numerator of h(x) under some assignment to x, it must also 
appear in the denominator of h(x), and vice versa. 

Lemma 4 (Square-Term Property): Given a coding vari- 
able x ee i, let /i(x) and / 2 (x) be the coefficients of x\ e , 
in TO a ;,(x)TO p9 (x) and m aq (x)m p i,(x) respectively. Then 
A(x) = / 2 (x). 

Proof: For any x ee >, define Qi = {(Pi, P 2 ) G P a bxP pg : 
x 2 ee , I f(Pi)t(P 2 )} and Q 2 = {(P 3 ,P 4 ) G V aq x V pb : x 2 ee , \ 
t(P 3 )t(P 4 )}. Consider a path pair (Pi,P 2 ) G Qi. Since the 
degree of x ee > in t(Pi) and t(P 2 ) is at most one, we must 



have x ee , | i(P x ) and x ee < \ t(P 2 ). Thus e, e' G Pi n P 2 . Let 
P/ , P 2 be the parts of Pi before e and after e' respectively. 
Similarly, define P 2 and P|. Then construct two new paths: 
P 3 = Pj U {e, e'} U Pi 2 and P 4 = P^ U {e, e'} U P| (see Fig. 
|3}. Clearly, t(Pi)t(P 2 ) - t{P 3 )t{P 4 ), and thus (P 3 , P 4 ) G Q 2 . 
The above method establishes a one-to-one mapping <j> : Qx — > 
Q 2 , such that for ^((Pi,P 2 )) = (P 3 ,P 4 ), t(Pi)t(P 2 ) = 
mXPt). Hence, /i(x) = ^ E^,*,)^ = 

^E(P 3 ,P 4)eS2 <(^)i(P4)="/ 2 (x). ■ 

V. Feasibility Condition of PBNA 

In this section, we provide the proofs of Theorems [T] [2] and 
3 (Main Theorem). 

A. 7y(x) Is Constant 

Proof of Theorem |7J As shown in |9|, when ?7(x) is 
constant and thus T = l2n+i> we set C = BA and Vi = 
(%)(2n+i)x(n+i)! where is an arbitrary variable, such that 
Eq. ((T|i is satisfied. Hence, if Pi(x) is not constant, 
and 33[ is satisfied. Thus R* is NA-feasible. Conversely, if 
Pi(x) is constant and P, = I 2n +i, SS\ is violated, and R* is 
not NA-feasible. ■ 

B. ?y(x) Is Not Constant 

Due to the importance of Vi, we first consider how to 
construct Vi which satisfies ([TJ. The construction of Vi 
involves solving a system of linear equations: 



r(z)(zC-BA) =0 



(7) 



where r(z) — (ri(z), ■ ■ ■ ,r n+ i{z)) € F 2 l ™ 1 (z). It is easy to 
see that zC — BA is a matrix on F 2 ™ (z). Assume ro(z) is a 
non-zero solution to Substitute z with ^(x 1 ), and we have 

ry(x 4 )r (ry(x l ))C = r (r ? (x i; ))BA 

Finally, construct the following precoding matrix 

Vf = (r^x 1 )) r^(x 2 )) ... tfWx*"- 1 ))) 

Apparently, Vi satisfies ([TJ. Hence, each non-zero solution to 
(j7]i corresponds to a row of Vi satisfying Q, Conversely, it 
is straightforward to see that each row of Vi satisfying ([T| 
corresponds to a solution to Q. 

Example 1: As an example, consider the case where n = 2 
and m — 2. Let a be the primitive element of F4 such that 



1 and a 



a - 



1 = 0. Moreover, let A = I 2 and 









^ 2 


a 


a 




- 


1 


1 


a 2 






V 1 


a 



Apparently, rank(C) = rank(B) = 2. It's easy to verify 
that r(z) = (a 2 z 2 + a, z + a, z 2 + az + a 2 ) satisfies 
equation (m). Thus, we substitute z with ry(x J ) and construct 
Vf = (r T ^(?7(x 1 )) r T (ry(x 2 )) •■■ r T (?7(x 5 ))). According 
to the above discussion, equation ([TJi is satisfied. ■ 
Using Q, we can derive the general form of Vj which 
satisfies Vi. 

Lemma 5: Any Vi satisfying ([TJ has the form Vi = 
GV^F, where VJ is defined in F is an (n + 1) x (n + 1) 



matrix, and G is a (2n + 1) x (2n + 1) diagonal matrix, with 
the (t, i) element being fi(i](x. 1 )), where fi(z) is a non-zero 
rational function in F 2 ™(z). Moreover, the (n + l)th row of 
FC and the 1st row of FBA are both zero vectors. 

Proof: See Appendix |B| ■ 
Lemma indicates that there is a direct relation between 
VI and the general form of Vi, which we use to prove that 
Eq. (|5j is also necessary for the feasibility of PBNA. 

Proof of Theorem BJ 1 The sufficiency of (|5]l was proved 
in J9J. Now assume p 4 (x) = g$$] & $n, where f(z) = 

J22 =0 a,kZ k and g(z) — Efe=o bkZ k - We will prove that for 
any Vi satisfying (jlj, ^ cannot be satisfied, thus R* is not 
NA-feasible. Apparently, if rank(Vi) < n+1, S%\ is violated. 
Thus, in the rest of this proof, we assume rank(Vi) = n + 1. 

By Lemma [5] Vi = GV*F, where F is an (n + 1) X 
(n+1) invertible matrix. The jth row of Vi is Tj = 
fj(r)(x j ))(l,r)(x j ),--- ,77 n (x^))F. Since the (n + l)th row 
of FC is zero, we have 

r j C = f j (r,(x j MlMx j ),--- ,»? n -V))H 
where H consists of the top n rows of FC and rank(H) = ra. 



Let a = (ao, ai, 



i T and b = (60, &i, 
F"^ and b' = FT 



1,2, we define a' 

r.a^/.-^x^a,^) 
= / J (r ? (x^))(l,r ? (x^") 

= pj( xJ ) r jCb' 



x b. 

,7fV))Fa' 
,77"(xi))a 



, 6„-i) J . For 
It follows 



" l (x»'))b 
- x (x*))Hb' 



Hence, the columns of (Vi PjViC) are linearly dependent, 
violating SS\. Similarly, we can prove the case of i — 3. ■ 
For the proof of the Main Theorem, we need to rearrange 

f(ri(x)) J I 

the ratio of rational functions ; ; k in Eq. (61 to a ratio of 
coprime polynomials with variables x. To this end, we use a 
property of polynomials stated in the following lemma. 

Lemma 6: Let F be a field, z is a variable and y = 
(yij V2, ■ ■ ■ ,Vk) is a vector of variables. Consider four non- 
zero polynomials f(z),g(z) £ ¥[z] and s(y),t(y) G F[y], 
such that gcd(/(z),5(z)) = 1 and gcd(s(y), i(y)) = 1. 
Denote d — max{d/, d g }. Define two polynomials in F[y]: 



and 



/3(y) = s(fgf)t d (y) 



Then 



gcd(a(y),/3(y)) = l. 

Proof: See Appendix [C] ■ 
The proof of the Main Theorem consists of two steps. In the 
first step, we use degree -counting technique and Linearization 
Property to reduce S n to the form { ""^"^^ }. In the second 
step, we use Linearization Property and Square Term Property 
to further reduce <S„ to the four rational functions in S'. 

Proof of the Main Theorem: Clearly, the necessity 
of (|6]l follows directly from Theorem [2] Now assume for 
i G {1,2,3}, pi(x) ^ S'. We will prove that p i (x) S n 
and thus R* is NA-feasible by Theorem |2] We only prove 
Pi(x) ^ S n . The other cases follow similar lines. By con- 



tradiction, assume there exists pi(x) 



44^ e S n . where 

/(z) = ^Lo a i z ' and = Si=0 & i^ such that a l b k ^ 
and gcd(/(«),ff(z)) = 1. Moreover, let p 1 (x) = ^ and 

^( x ) = f{$> where g cd ( M ( x )> w ( x )) = gcd(s(x),t(x)) = 1. 
Let d = max{fc, I}. Define the following two polynomials 
a(x) = f(r](-x))t d (-x) and /3(x) = g(^(x))i d (x). According 
to Lemma [6] gcd(a(x), /3(x)) = 1. Thus, we have a(x) = 
cit(x), and /3(x) = cv(x), where c <E and c 7^ 0. 

According to Lemma [3] there exists an assignment to 
x under which it(x) and v(x) are transformed into either 

U(x ee >) = CiX ee ' + C , v(x ee >) = C 2 Or u(x ee >) = C 2 , 

^(a^ee') = Cix ee > + c . We only consider the first case. 
The proof for the other case is similar. In this case, a(x) 
and /3(x) are transformed into a(x ee ') = ccix ee > + ceo and 
f3(x ee ') = CC2 respectively. 

First, we prove that both t(x ee i) and s(x ee /) are non-zeros. 
Assume t{x ee i) = 0. If k ^ I, at least one of a(x ee >) and 
j3(x ee ') equals zero, which is impossible. On the other hand, 
if k = I, we have a(x ee i) — akS k (x ee ') and f3(x ee >) = 
bkS k (x ee >). It follows that cc\x ee i +ccq = akb~ k 1 cc 2 , which is 
impossible. Thus we have proved that t(x ee i) 7^ 0. Similarly, 
we can also prove that s(x ee >) 7^ 0. 

We then prove that d = 1. By contradiction, assume d > 2. 
We first consider the case where I < k and thus d = k. In this 
case, we have 

a(x ee >) = ' ^a j t k ~ 3 {x eel )s ] {x ee >) = ccix ee , + cc 

(3(x ee >) = _ Q bjt k ~° (Xee^S 1 {x ee >) = CC 2 



Define qi(x) 



Pl(x) 

consider the following cases. 



mn (x)m 2 3(x) 
mi 3 (x)m 2 l(x) ' 



For d 



1, we 



Case I: 



b +b 1 z 



where aiao&i&o 7^ 0, and ao&i 7^ 

ai b . For this case, we have P i(x ee >) = bo+biPufc'^uX"') ' 
It immediately follows 

a C2 - 60C0C2 

gi (a; ee ' ) = 



boCxC2X e 



b\c\x 2 ee , - axc^Xee* + 6iCq - aic c 2 

Denote wi(a; e e') = aoc 2 ~~ fy)CoC2 — b$c\C2X ee i and vi^ee') = 
b\c\x 2 ee , — a\C\C2X ee ' + &iCg — aiCoC2. Assume wi(iEee') | 
v i( x ee') an d thus a; ee ' = aoC fc "^ oCo i s a solution to 1>i (iEee' ) = 

^^■(ao^i + ai6g) 7^ 0. Hence, 



0. However, vU a ° c ?- b ° c ° 



ui(x ee ')\ vi(x ee >). Thus, by the definition of <?i(x) and 
Lemma [4 s^ e , must appear in w^^ee')' which contradicts 
the formu ation of u\(x ee i). 

Case II: 44 = a °+° lZ , where a ai6 7^ 0. Similai- to Case 
I, we can derive 



qi (x e 



a c\ 



b i c i x te' - aicic 2 x e 



he? 



5 - aic c 2 



which contradicts Lemma |4] 



Case III: 



ZM _ 



b +b 1 z 



where aibobi 7^ 0. Thus 



a! 6 "ii 3 (x)m 2 i(x) 
61 61 mn(x)m 2 3(x) 
Since the coefficient of each monomial in m n (x)m23(x) and 



Case V: 4^ 



Pi(x) £ 5'. 

Assume s(x ee >) = ^=0 s j x le< and K x ee>) = Ej=o^ x ee" Case IV: 
where s r t r / 7^ 0. Thus r = d s and r' = d t and max{r, r'} > 
1. Note that the degree of a; ee ' in t k ~i {x ee ')s 3 (x ee ') is kr' + 
j(r — 7-'). We consider the following two cases: 

Case I: r ^ r'. If r > r', d a = kr > 2, contradicting that 
d a = 1. Now assume r < r' . Let and I2 be the minimum 
exponents of z in /(z) and g(z) respectively. It follows that 
d a = kr' - h{r' - r) = 1 and dp = kr' - l 2 (r' - r) = 0. 
Clearly, l 2 > due to dp = 0. If r > 0, /fcr' - Z 2 (r' - r) > 
kr' ~l 2 r' > 0, contradicting = 0. Hence, r = 0, and I2 = k 
due to d^j = 0. Meanwhile, d Q = (k—li)r' = 1, which implies 
that ?i = k — 1 and r' = 1. Thus, z fe_1 is a common divisor 
of f(z) and g(z), contradicting gcd(/(z), g(z)) = 1. 

Case II: r = r'. Since d a = 1 and dp(x ee >) = 0, all the 
terms in a(x ee > ) and [3(x ee > ) containing a;^/ must be cancelled 
out, implying that 



mi3(x)rri2i(x) equals one, it directly follows ^ = — ^ = 
5^ = 1. This indicates that pi(x 



61 bi 
contradicting 



/M _ 

qi {Xee>) 



, where a b bi 7^ 0. It follows that 



g c| - 6 c C2 - 6 cic 2 x e 

2™2 



&1C5 + 6icja; 



Similar to Case I, this also contradicts Lemma [4] 

9 f, where clq 7^ 0. Hence, qi(x ee <) 



2 2 x-g , contradicting Lemma 



Case VI:g} = ao- 



Pl(x) = a + ai 



aiz, where a^ai 7^ 0. Thus, it follows 

m 3 i(x)mi2(x)m 2 3(x) 



Similar to Case EH, a\ 



Case VII: 



IM _ 
9(2) 



m 2 i(x)TO3 2 (x)mi 3 (x) 
ap = 1, contradicting pi(x) ^ 5'. 



j=o j=o 



( *' 



^ ' r 











Hence z — is a common divisor of /(z) and ,g(z), contra- 
dicting gcd(/(z),5(z)) = 1. 

Therefore, we have proved d = 1 when Z < fc. Using similar 
technique, we can prove that d = 1 when I > k. 



a±z, where ai 7^ 0. Similar to Case III, 
Pi(x) = ry(x), contradicting pi(x) ^ iS'. 

In summary, we have proved that if pi (x) ^ S', p L (x) ^ S n 
and hence R* is NA-feasible by Theorem [2] ■ 

C. Some mij(x) =0 (i ^ j) 

In this case, since the number of interference terms is 
reduced, at least one of sd\ , s^2 , ^\ is removed, and thus the 
restriction on Vi imposed by Eq. ([TJ vanishes. Therefore, 
we can choose Vi freely, and the feasibility condition of 
PBNA is greatly simplified. For example, assume to 2 3( x ) = 
and all other transfer functions are non-zeros. Hence £#2 is 
removed. Meanwhile, SS^, 3§' 2 > ^3 remain the same. Similarly 



to Theorem[T] we can set Vi = (0y)(2n+i)x(n+i)> where 9^ 
is an arbitrary variable. It is easy to see that R* is NA-feasible 
if and only if Pi(x) is not constant for every i E {1,2,3}. 
Using similar arguments, we can discuss other cases. 

VI. Checking the Feasibility of PBNA 

For a given graph, checking the feasibility of PBNA is now 
reduced to checking whether p;(x) does not have any of the 
four simple functional forms in Eq. for i G {1, 2, 3}. This 
is a multivariate polynomial identity testing problem. In gen- 
eral, polynomial identity testing is a Randomized Polynomial 
time (RP) problem and is currently not known to be in P. 

To check whether PiCx) ^ 1, we use Ford-Fulkerson Algo- 
rithm, as per Lemma [2] To check whether p;(x) ^ ??(x), we 
define gi(x) = ^ ;vjjj and consider g;(x) ^ 1. Therefore, Ford- 
Fulkerson Algorithm can be used to check this condition as 
well. For the last two conditions (Pi(x) ^ {l + ?;(x), ^^-fi })' 
it is still not clear what is their interpretation in terms of 
graph structure. A counter example is shown in Fig. |TJb). 
Nevertheless, we can still check the conditions by evaluating 
the rational functions through T random tests: 

for k = 1 to T do 

Assign random values to x, denoted by Xo 

If Pi(x ) f {1 + r/(x ), „^+], }i return success 

end for 



Return failure (i.e., is violated) 



Let L denote the maximum distance from any sender to any 
receiver in the network. Using Lemma 4 of p), we can upper- 
bound the probability of error as follows. We consider the case 
of i = 1. Other cases follows along similar lines. Note that 
Pi(x) = 1 + rj(x) is equivalent to the following equation: 

/(x) = m 3 i(x)TOi2(x)TO2i(x)mi 3 (x) + 
mn (x)mi2 (x)m 32 (x)m 13 (x)+ 
mii(x)m 31 (x)m 12 (x)m 23 (x) = 

Since the maximum degree of any variable x ee ' in a transfer 
function is at most one, the total degree of each term in /(x) 
is at most 4L. For each random test, the probability of error 
in checking if pi(x) ^ 1 + 77 (x), denoted by 81, can be upper 
bounded by using Lemma 4 of |3J: Si = Pr(/(x ) = | 
/(x) 7^ 0) < 1- (l - 57^2) L ■ Hence, the total probability of 
error in checking if pi(x) ^ 1 + ?7(x) is Pi(Error) = Sj < 
[1 — (1 — 2 „}-2 ) L ] T ■ Similarly, the total probability of error in 
checking if pi(x) ^ v ?x)+i < denoted by P2(Error), is upper 



bounded by P2(Error) < [1 — (1 



\L1T 



. In summary, the 



total probability of error is upper bounded by: 

P(Error) < Pi(Error) + P2(Error) 



< 2 



1-1- 



1 



Thus, the error can be made arbitrarily small for suffi- 
ciently large m and T. The running time of the algorithm 
is 0(T\E\Di n ), where Di n is the maximum in-degree of any 
node in the network. 



VII. Conclusion 

In this paper, we study the feasibility of PBNA for three 
unicast sessions. We first prove that the set of conditions 
proposed by |9| are also necessary for the feasibility of 
PBNA with respect to any valid precoding matrix. Then, we 
reduce this set of conditions to just four conditions, using two 
graph-related properties along with a simple degree-counting 
technique. This reduction enables an efficient algorithm for 
checking the feasibility of PBNA. 

Appendix A 
Proofs of Graph Properties 

The following lemma is used in the proof of Lemma [2] 
Lemma 7: Let (Pi,P2) E Tab x 'Ppq- Then, there exists 
(P 3 ,P 4 ) € V aq x V pb such that t(Pi)t(P 2 ) = t(P 3 )t(Pi) if 
and only if Pi n P 2 + 0. 

Proof: First, Assume Pi n P% 7^ 0. Pick an arbitrary edge 
e € Pi PI P2. Let P\ and Pf be the path segments along Pi 
before and after e respectively. Similarly, we can define P\ 
and P|. Construct P 3 = P\ U {e} U P x 2 and P 4 = P\ U {e} U 
P|. Hence, it is easy to see that (P 3) P4) 6 V aq x V v b and 
t(Pi)t(P 2 ) = t(P 3 )t(P 4 ). 

Now assume Pi n P2 = 0. By contradiction, assume 
there exists (P 3 ,P 4 ) G V aq x V pb such that t{Pi)t{P 2 ) = 
t(P 3 )t(P 4 ). Clearly, P x U P 2 = P 3 U P 4 . Then, there exist 
e,e' G P 4 such that head(e) — tail(e') and e £ Pi, e' E P2. 
Hence, x ee ' \ £(P 3 )i(P 4 ) but x ee > \ i(Pi)i(p2), contradicting 
our assumption. ■ 
Proof of Lemma |2| Assume m a {,(x)rn pg (x) 7^ 
m ag (x)m p fc(x). Thus there exists (Pi,P2) E V a b x V vq such 
that for any (P 3 ,P 4 ) E V aq x V pb , t(Pi)t(P 2 ) ^ t(P 3 )t(P A ), 
or vice versa. By Lemma [7] P 1 n P 2 = (P 3 n P 4 = 
for the other case). On the other hand, if there exists dis- 
joint path pair (P^Pj) E V a b X P pg , t(Pi)t{P 2 ) is absent 
from TO a9 (x)m p b(x). Moreover, there is only one term in 
m a b(x)m pg (x) which equals i(Pi)£(p2). Thus £(Pi)£(P 2 ) 
doesn't vanish from m ab {x)m pq {x). Hence m a b(x)m pg (x) 7^ 
m ag (x)m p b(x). Similarly, the theorem holds for the other 



Appendix B 
General Form of V x 

The following lemma shows that given any full-rank ma- 
trices A, B and C as defined in srf{, ^ > we can always 
find a non-zero solution to dTli, and thus construct a precoding 
matrix Bi which satisfies Q. 

Lemma 8: Equation ^ has a non-zero solution in F 2 l i 1 [z] 
in the form of r(z) = (l,z,z 2 ,--- ,z n )F, where F is an 



(n + 1) X (n + 1) matrix in ¥2^. Moreover, any solution to 
(|7| is linearly dependent on (1, z, ■ ■ ■ , z n )F. 

Proof: Denote D = BA. First, we will prove that 
rank(zC — D) = n. Let Cj and denote the ith col- 
umn of C and D respectively. Hence, , • • • , c„ are lin- 
early independent and so are di, • • • ,d ra . Assume there ex- 
ist fi(z),--- ,f n (z) E ¥ 2 m(z) such that Yh=i fi( z )( zc i 
dj) = 0. Without loss of generality, assume fi(z) 



h(z) 



for i E {1,2,- •• , n}, where gi(z),h{z) E F2m[z]. Thus, 



Yn=i9i{ z ){ zc i ~ d i) = 0. Let k = max ie {i )2) ... , n }{d gi } 
and assume gi(z) — X)f=o a ' i2 '- Then, it follows 



^2g l (z)(zc l - dj) = ^^(a^ i z /+1 c i - aj^'d;) 

(=0 i=l 
n k — 1 n 

1 ^2 a>k,iCi + ^2 z ' +1 ^2( a ^i c i ~ a i+id d i) 



i=X 



1=0 i=l 



- ^2 a o,t d » = 

i=l 

Therefore, the following equations must hold: 

n n 

^2 ak ' iCi = ^2 a °' ,:d,; = 

i=l i=l 
n 

J2( a u c * - a i+i,i d i) = w g {o, 1} 

i=l 

Thus a;,i = for any i £ {1, ••• ,n},l £ {0, ■ ■ • , k}, 
implying fi(z) — 0. Hence, rank(zC — D) = n. 

Then, there must be an n x n invertible submatrix in zC— D. 
Without loss of generality, assume this submatrix consists 
of the top n rows of zC — D and denote this submatrix 
by E„ + i. Let b denote the (n + l)th row of zC — D. In 
order to get a non-zero solution to equation Q, we first fix 
r n+ i(z) = — 1. Therefore, equation (j7| is transformed into 
(ri(z),--- ,r n (z))E„+i = b. For i £ {1,2,- •• ,n}, let E, 
denote the submatrix acquired by replacing the ith row of 
E„ + i with b. Hence, we get a non-zero solution to (|7]): 

detEx detE ra 
v ; v dctE„+i' ' detE n+ i' ' 

Moreover, f(z) = (det E x , ■ ■ • , det E n , — det E n+1 ) is also a 
solution. Also note that the degree of z in each det E; (i £ 
{1, 2, • ■ • , n + 1}) is at most n. Thus, r(z) can be formulated 
as (1, z, z 2 , ■ ■ ■ , z n )F, where F is an (n+ 1) x (n + 1) matrix 
in F2™. Since rank(zC — D) — n, all the solutions to equation 
(|7]i form a one-dimensional linear space. Thus, all solutions 
must be linearly dependent on f (z). ■ 
Proof of Lemma^ Let be the ith row of Vi, which 
satisfies equation ([T]). According to Lemma|8] r, must have the 
form /i(?7(x I ))(l, ?y(x J ), • • • ,r; n (x l ))F, where /i(z) is a non- 
zero rational function in ¥2^(z). Hence, Vi can be written 
as GV^F. 

According to Lemma [8] equation (|7]i can be rewritten as 
follows: 

(z, z 2 , ■ ■ ■ , z" +1 )FC = (1, z, ■ ■ ■ , z")FBA 

The right side of the above equation contains no z n+1 , and 
thus the (n + l)th row of FC must be zero. Similarly, there 
is no constant term on the left side of the above equation, 
implying that the 1st row of FBA is zero. ■ 

Appendix C 
Results on Multivariate Polynomial 



i.e., the vector consisting of all variables in y other than j/j. 
Note that any polynomial /(y) G F[y] can be formulated as 

/(y) = /o(y») + hfa)yi h 1- f P fa)Vi 

where ft fa) G F[y 2 ] for j G {0, 1,- •• ,p} and / p (y,-) ^ 0. 
Let F(yj) denote the field consisting of all rational functions in 
the form of ^4, where ufa),v(yi) G F[yj]. Because F[y.;] 
is a subset of F(yj), /(y) can also be viewed as a univariate 
polynomial in the ring F(y^) [j/j]. For any h(y) £ F[y], we 
use h(yi) to denote its equivalent counterpart in F(yj)[^]. 
To differentiate these two concepts, we reserve the notations, 
such as "|", "gcd" and "lcm'j^] for field F, and append "1" as a 
subscript to these notations to suggest they are specific to field 
F(yj). For example, for f(y),g(y) G F[y] and u{yi),v(yi) £ 
r (yi)[Vil g(y) I /(y) means that there exists h(y) £ F[y] 
such that f(y) = h(y)g(y), and u{yi) |i v(y l ) means that 
there exists w(yi) £ W\yi\(yi) such that v(yi) = w(yi)u(yi). 
Similarly, gcd(/(y), g(y)) is the greatest common divisor of 
f(y) and g(y) within F[y], and gcA 1 (u{y l ) 1 v{y i )) is the 
greatest common divisor of u(yi) and v(yi) in F(y;) [?/;]. 

In general, each polynomial h(yi) £ F(yj)[?/j] is of the 
following form 



a*ofa) , ai(yi) 



-Vi 



q P (y») P 
1 y% 



&o(yi) 61 (yi) ^(yO' 

In the above formula, for any j G {0, 1, • • • ,p}, 
aj(yi),bj(yi) £ F[yj], ^(y,) 7^ 0, gcd(a i (y i ),6 J -(y i )) = 1, 
and a p (yj) 7^ 0. Note that for any y\ which is absent in 
h(yi), we let a,-(yi) = and &j(yi) = 1. Define the following 
polynomial 

Wj(yi) = lcm(6o(yi),6i(y<), • • • ,b p (yi)) 

Thus, HhiYi) G F[yi] and n h (yi)h(yi) £ F[y]. 

Lemma 9: Assume g(yj) G F[y.j] and /(y) G F[y] is of 
the form /(y) = ^=0 /j(yi)»J. where /j(yi) e F[y,]. 
Then g(y,) | /(y) if and only if g(yj) | ^(yj for any 
jG{0,l,--- ,p}. 

Proof: Apparently, if g(y{) \ ft fa) for any j £ 
{0, 1, • • • ,p}, gfa) I /(y). Now assume ^(yi) | /(y). Thus 
there exists /i(y) G F[y] such that f(y) = g(yi)h(y). Let 
My) = Y, P j=o h j(yi)yi' Hence, it follows that ft fa) = 
hjfa)gfa) and thus gfa) \ ft fa). ■ 
The following result follows immediately from Lemma [9] 
Corollary 1: Let <?(y,) and /(y) be defined as Lemma [9] 
Then gcd(g(y,), f(y)) = gcd(gfa),f fa), ■■■ , f p fa)). 

Proof: Note that any divisor of g(yi) must be a poly- 
nomial in F[yi]. Let dfa) = gcd(gfa), /(yl) and d'(y,) = 
gcd(g(y 4 ),/o(yi)," - ,f P fa))- By Lemma[9| rf(y 4 ) | ft fa) 
for any j G {0, 1, • • ■ ,p}, implying that dfa) \ d'fa). On 
the other hand, d'fa) \ f(y), and thus rf'(yi) | d(yj). Hence, 
dfa) = d'fa). ' ' m 

Corollary 2: For t £ {1,2, ••• ,s}, let / t (y) G F[y] be 
defined as f t (y) = J2% ftjfa)Vi' where ftjfa) G F[y,]. 



Let y = (yi, 2/2) " " " j Vk) be a vector of variables. For any 
i £ {1,2, ••• ,k}, define y 4 = (yi,--- ,yi-%, yi+i, ■■ ■ ,yk), 



5 We use lcm(/(2;), g(x)) to denote the least common multiple of two 
polynomials f(x) and g(x). 



Let g(yi) <= F[y j; ]. It follows 

gcd(ff(y i ),/i(y), • • • ,/ t (y)) 

=g cd (5(y l ), /io(yi), • • • )/i P1 (yi),--- , 

/ao(y»)) ■ • • i/«P,(yt)) 
Proof: We have the following equations 

gcd(g(yi),/i(y) } --- ,/ t (y)) 
=gcd(sf(yi), /i(y), • • ■ ,s(yi), / t (y)) 
=gcd(gcd( fl (y i ), /i(y)), • • • , gcd( 5 (y 4 ), / a (y))) 
=gcd(5(yi),/io(yi),--- ,/i Pl (yi),--- , 

ff(yi),/ao(yi),--- ./«p,(yi)) 
=gcd(5(yt)>/io(y*))--- ,/i P i(yi),--- > 
/ s o(yi),-- - ./ap.(yi)) 

■ 

Lemma 10: For £ € {1,2, ••• ,s}, let a t (y),b t (y) € F[y] 
such that &t(y) 7^ and gcd(a t (y), 6 t (y)) = 1. For £ e 
{1,2,- •• ,s}, let «t(y) = lcm(6i(y),--- ,6 t (y)). Then we 
have 

gcd ai(y)— ,a s (y)—— ,v s (y) 



'My)' '" BW/ 6 s (y)'"' -■' ' ' 

Proof: We use induction on s to prove this lemma. Appar- 
ently, the lemma holds for s = 1 due to gcd(ai(y), 61 (y)) = 
1. Assume it holds for s — 1. Thus it follows 



gcdlai(y) 
=gcd^ ai (y) 
=gcd( ai(y) 



= gcd^ai(y) 
= gcd( ai(y) 



gcd ai(y) 



^(y) 

61 (y)'' 

^(y) 

61 (y)'' 

My) 

My)'' 

^(y) 

61 (y)'' 

My) 

My)'' 

My) 



/ N^(y) / \ 
My) 

,My)^4,My)^ (y) 



My) 

,gcd(a s (y),6 s (y)) 



My) 
My) 



,a«-i(y) 
,a s -i(y) 



b s (y) 
v s {y) v s (y) 



cd v s -i(y) 



Mi(y)' My) 
My) 
6-1 (y)' 
My)' 



My) 



=gcd 



My) 
My) 



, \ My) , \ My) 

> a «- 1 (y) 1 — ' (y) 



gcd 



w«-i(y) 
My) 



gcd ai(y) 



w s -i(y), 



Mi(y) 
^s-i(y) 
My) 
My) 



Os-i(y) 



My) 
««-i(y) 
6«-i(y) 



My) 



u«-i(y) 



,w«-i(y)i 



(<0 



My) 



gCd Vgcd(u s _i(y),6 ;j (y)) 
gcd(l,«._ 1 (y)) = l 



My) 
My). 

v*-i(y)) 



« s -i(y) 



jcd(w s _i(y),6 s (y)) 



In the above equations, (a) is due to gcd(a s (y), b s (y)) = 
1; (b) follows from the fact that 

thus ^g} = gcd(^_ 1 (y) 



1 Vs-i(y) and 
b / >' ); (c) follows from the 



inductive assumption; (d) is due to the equality: v s (y) 
lcm(« s _ 1 (y),6 s (y))- 



gcd(u s _i(y),fc s (y)) ■ 

Corollary 3: For j G {1, 2, • • • , s}, let £(3/,-) € F(y*)[sfc]. 
Define u(y ( ) = lcm(/x /l (y* )_,-■■, /i/, (y,)) and fj(y) = 
v(yi)fj(Vi). Thus gcd(u(yO,/i(y), ■ ■ ■ ,/ s (y)) = 1 
Proof: Assume fj(yi) has the following form: 



,(,_ My*) , My*),. 
&io(y») »ii(y») 



a 3Pj(yi) Pj 

KM) 



where for any j € {1,2, ••• ,s} and £ € {0,1,- •• ,f>j}, 
a jt(yi),bj t (yi) e F[yJ, b,- t (yi) 7^ and 
gcd(a Jt (y i ),6j t (y i )) = 1. Apparently, v(y t ) is 
the least common multiple of all M(y$)' s - Define 



Ujt(yi) 



My.) 



6 F[yJ. Hence, we have 



/i(y) = Etio a J*(yi) u jt(y<)2/*' Then k follows 

gcd(«( yi ),/i(y)),--- ,/ s (y)) 

= gcd(w(y l ),aio(y 4 )wio(y l ), • • • , ai Pl ( yi )ui pi (y*), • • • , 
a S Q(yi)u s o(yi), ■ ■ ■ ,a S p s (yi)us Ps (yi)) 

where (a) is due to Corollary [2] and (b) follows from Lemma 

M ■ 

Generally, the definitions of division in F[y] and F(yj)[t/i] 
are different. However, the following theorem reveals the two 
definitions are closely related. 

Theorem 4: Consider two polynomials /(y),y(y) G F[y], 
where g(y) ^ 0. Then g(y) \ f(y) if and only if gfa) | x /(j/j) 
for every i € {1, 2, • • • , fc}. 

Proof: The division equation between /(j/j) and is 
as follows 



f(Vi) = h i {y i )g{y l ) +n{yi) 



(8) 



where h i (y i ),r l (y i ) € F(y^)[^], and either n(yi) = or 
d ri < d g . Due to the uniqueness of Equation (pjj, /(y) | 5(y) 
immediately implies that for any i € {1, 2, • • • 7^}, ^(j/i) = 
and thus |i /(?/,). 

Conversely, assume for every i 6 {1, • • • , fc}, <7(j/j) |i /(yi) 
and hence r»(j/i) = 0. Denote ^(y) = (ihi(yi)fk(yi). Clearly, 
' l i(y) € F[y]. Then, the following equation holds 

Mfci(y»)/(y) = h l {y)g{y) 

By Corollary [3] gcd^fy,-), /i,(y)) = 1. Thus, ^(y^) | 
5 (y). Define g(y) = ^g^. By Lemma g g(y) E F[y]. 

Define w(y) 



^ € F[y]. It follows that 



u(y) = 



gcd(/(y), s (y)) 



gcd(/(y), ff (y)) 

^(y»)g(y) 

gcd(/ii(y)5(y),Mhi(yi)fl(y)) 
Mfe,_(yQg(y) 

ff(y)gcd(/i i (y),/i /li (y i )) 

Mfe«(yi)g(y) 

My) 
Vhi (y») 



Note that variable t/j is absent in u(y). Because x/i can be 
any arbitrary variable in y, it immediately follows that all the 
variables in y must be absent in u(y), implying that u(y) is 
a constant in F. Hence g(y) | /(y). ■ 

Theorem 5: Let /(y),g(y) be two non-zero polynomi- 
als in F[y]. Then gcd(/(y), g(y)) = 1 if and only if 
gcd 1 (/(y J ), g{vi)) = 1 for any i G {1,2,- • • ,k}. 

Proof: First, assume for any i G {1,2, ••• , k}, 
gcd 1 (/(j/i), g{yi)) = 1. We use contradiction to prove that 
gcd(/(y),s(y)) = 1. Assume u(y) = gcd(/(y), g(y)) is 
not constant. Let y^ be a variable which is present in u(y). 
By Theorem [5] u(yj) |i /(y,) and u(j/i) |i which 
contradicts that gcA 1 (f{y i ),g{y i )) = 1. 

Then, assume gcd(/(y), fl'(y)) = 1. We also use 
contradiction to prove that for any i G {1,2, ••• , k}, 
gcd 1 (f(y i ),g(y i )) = 1. Assume there exists i G {!,-■■ ,k} 
such that v(yi) = gcd 1 (/(yj), g(yi)) is non-trivial. Define 
w(y) = fJ-v(yi)v(Vi) G F[y]. Clearly, |i /(y*) and 

w(yi) |i Thus, there exists p{yi),q{yi) G F(y i )[y i ] such 
that 

/(s/i) = w(y l )p(y i ) g{yi) = w(y % )q{yi) 

Let sjyi) = lcmO^yi),/^^)). Define p(y) = s(yi)p{y,) 
and q(y) = s(y l )q(y i ). Apparently, p(y),q(y) G F[y]. It 
follows that 

s (y»)/(y) = w(y)p(y) s(yi)g(y) = w(y)q(y) 
Then the following equation holds 

s(y<)gcd(/(y),5(y)) = w(y)gcd(p(y),g(y)) 



Due to Corollary [3] gcd(s(yi),; 
gcd(s(yi),_p(y),q(y)) = 1. Hence 



w(y) 



"(y) 



;cd(p(y),g(y))) 
sf-Vi) I w(y). 



According to Lemma 



Let 



u>(y) is a non- 
?cd(/(y),y(y)), 



trivial polynomial in F[y]. Thus, w(y) 
contradicting gcd(/(y), g(y)) = 1. ■ 
Lemma 11: Consider two non-zero polynomials in F[z], 



f(z) = aa+aiz + - ■ - + a p z p and g{z) = bo + biz + - ■ 
where cii,bj G F for i G {0, 1, • • • ,p},j G {0,1. 
a p b q ^ 0, p > q and gcd(/(z), g{z)) — 1. Let s(x),i(x) 
two non-zero polynomials in ¥[x] such that gcd(s(x),t(x)) 
1. Define the following polynomials in F[x]: 



,q}, 



p{x) 



f 



9 



t(x) 
s(x) 



a k tP- k (x)s k (x) 
b k t p - k (x)s k (x) 



t(x) 

Then gcd(a(x),/?(x)) = 1. 

Proof: Assume it)(a;) = gcd(a(x), /3(x)) is non-trivial. 
Thus we can find an extension field F of F such that there 
exists xo G F which satisfies w(xo) = and hence a(xo) = 
f3(xo) — 0. In the rest of this proof, we restrict our discussion 
in F. Note that gcd(/(z), g{z)) = 1 and gcd{s(x),t(xj) = 1 
also hold for F. Assume t(xo) = and thus x — xo | t(x). 
Since gcd(s(x),t(x)) — 1, it follows that x — x n \ s(x) and 
thus s(xq) 0. Hence, a(xo) = a p s p (xo) ^ 0, contradicting 
that a(xo) = 0. Hence, we have proved that t(xo) ^ 0. Then 



we have 

's(xo)A 
t(x )J 



f 



a(x ) 
tP(x ) 
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S(X ) \ 

t(x )J 



tP(x ) 



which implies that z 



t(xo 



is a common divisor of f(z) and 



g(z), contradicting gcd(/(z), g(z)) = 1. ■ 
Proof of Lemma [6[ Note that if we substitute F with 
F(yi) and gcd with gcdj^ in Lemma 11 the lemma also 
holds. Apparently, f(z),g(z) G F(yjj[z]. We will prove 
that gcd 1 (/(z), g(z)) = 1. By contradiction, assume r(z) = 
gcd 1 (/(z), 5 (z)) G F( yi )[z] is non-trivial. Let f(z) = 

and g(z) = Clearly, f(z) and g{z) are both non-zero 

polynomials in F(yj)[z]. Then we can find an assignment to 
y,:, denoted by y*, such that the coefficients of the maximum 
powers of z in r(z),f(z) and g(z) are all non-zeros. Let 
f{z) denote the univariate polynomial acquired by assigning 
yi = y* to r(z). Clearly, r(z) is a common divisor of f(z) and 
g(z) in F[z], contradicting gcd(/(z), g(z)) = 1. Moreover, due 
to gcd(s(y),t(y)) = 1 and TheoremBJ gcd 1 (s(y i ), ifa)) = 1. 
Thus, by Lemma [IT] gcd 1 (a(yi),/3(yi)) = 1. Since i can be 
any integer in {1, 2, • • • , k}, it follows that gcd(a(y), /3(y)) 
1 by Theorem |5] ■ 
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